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^ ABSTRACT 
D 

^ ■ In the supersymmetry (SUSY) models, because of the mass difference 

between the left- and right-top squarks, the supersymmetric QCD in- 
^ \ teractions can generate parity violating effects in the production of tt 

pairs. We show that SUSY QCD radiative corrections to the parity vi- 
olating asymmetry in the production rates of the left- and right-handed 
top quarks via the qq tt process can reach about 3% at the Fermilab 
Tevatron with ^J~s = 2 TeV. This could be observable with an integrated 
luminosity larger than 2fb _1 . We also show that in these models the 
SUSY QCD radiative corrections to the ti production rate are small 
unless gluinos are very light, of the order of 1 GeV. 



1 Introduction 



In a recent paper we studied the parity violating asymmetry induced from the 
super symmetric electroweak (SUSY EW) and Yukawa (SUSY Yukawa) corrections at the 
one loop level. Two classes of supersymmetry (SUSY) models were considered: the mini- 
mal supergravity (mSUGRA) models 0] and the minimal supersymmetric models (MSSM) 
with scenarios motivated by current data || Q . After sampling a range of values of SUSY 
parameters in the region that might give large contributions to the parity-violating asym- 
metry A, and which are also consistent with either of the above two classes of models, 
we found that the asymmetry A due to the one-loop SUSY EW (amj/m^) and SUSY 
Yukawa corrections for the production process qq — > g — > ti at the upgraded Tevatron 
is generally small, less than a few percent. However, the sign can be either positive or 
negative depending on the values of the SUSY parameters. (The effect from the Standard 
Model (SM) weak corrections to this asymmetry is typically less than a fraction of percent 

!!•) 

In the supersymmetric Standard Model, some superparticles experience not only the 
electroweak interaction but also the strong interaction. Although the SM QCD interaction 
respects the discrete symmetries of charge conjugation (C) and parity (P), the SUSY QCD 
interactions for superparticles, in their mass eigenstates, need not be C and P invariant. 
(Needless to say, in the strong interaction eigenstates, the SUSY QCD interaction is C- 
and P-invariant.) For either the mSUGRA or the MSSM models, the masses of the left- 
stop (the supersymmetric partner of the left-handed top quark) and the right-stop can be 
noticeably different due to the large mass of the top quark. This is a general feature of 
the supersymmetry models in which the electroweak symmetry is broken spontaneously 
via radiative corrections. Since both the left-stop and the right-stop contribute to the 
loop corrections for the ti pair production process qq, gg — > ti, the different masses of the 
top-squarks will induce a parity violating asymmetry. It is this effect that we shall study in 
this paper. Because the ti pairs are produced predominantly via the QCD process qq — > ti 
at the Tevatron (a pp collider with CM energy a/s = 2TeV), we shall concentrate on 
the SUSY QCD corrections for the qq fusion process. We show that the parity violating 
asymmetry at the Tevatron induced by SUSY EW, SUSY Yukawa and SUSY QCD 
effects could add up to more than 3% for some of the SUSY models. The SM tree level 
cross section is 4.34 pb for a 176 GeV top quark. The SM QCD correction increases the rate 
by about 60% |7J, while the SM EW correction is less than a percent ||. It is expected 
that there are about 2300 fully reconstructed 6-tagged ti events in the lepton plus jets 
mode collected by the two experimental groups with a 2 fb -1 integrated luminosity |J. 
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This amounts to a sie nal at ~ 90% c.l. (confidence level) with 2 fb" 1 , or 99% c.l. with 10 
fb~ . Thus, a study of A at the Tevatron could yield information about the allowed range 
of SUSY model parameter space. 

2 SUSY QCD Corrections and Parity Violation 



I. Squark mixings 

In the MSSM the mass eigenstates qi and q 2 of the squarks are related to the (strong) 
current eigenstates q^ and q R via the mixing angle 6q by 

qi = q L cos 6q + q R sin 6<j, q 2 = -q~L sin 6g + q R cos 9g. (1) 

For the top squarks, the mixing angle % and the masses m~ tl 2 can be calculated by diago- 
nalizing the following mass matrix ||, 

M 2 _ ( M l m ^LR \ 

1 \m t m LR M? r J' 



M h = m h + m t + (3 ~ 3 sin2 9w "> cos ( 2 ^) m i 
2 

M? r = m~ R + m\ + - sin 2 ^ cos(2/5)m|, 

= -/icot/3 + A t , (2) 

where m 2 and m 2 are the soft SUSY breaking mass terms of the left-stop and the 
right-stop, n is the coefficient of the H1H2 mixing term in the superpotential, A t is the 
parameter describing the strength of soft SUSY breaking trilinear scalar interaction tx,t^if 2 , 
and tan/5 = v 2 /vi is the ratio of the vacuum expectation values of the two Higgs doublets. 
9w is the weak mixing angle, and mz is the mass of the Z boson. 

II. Renormalized amplitudes and the asymmetry 

The effects of parity nonconservation can appear as an asymmetry in the invariant mass 
(M t t) distributions as well as in the integrated cross sections (ex) for ti and t R production. 
The integrated asymmetry, after integrating over a range of M t t, is defined by IQ 

A = ^^ = ^^, 0) 
N R + N L ctr + gl 

where iV^ and Nl are the numbers of right-handed and left-handed top quarks. Thus, 
&r = ctrl + vrr (ol = &LR + o"ll) is the cross section for producing a ti pair in which the 
top quark is right- (left-) handed, and the top antiquark is either left- or right-handed. 
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Some of the one loop scattering amplitudes of qq — > it were already presented in 



Refs. [|I0| , |TT|j for calculating the total production rates of ti pairs. To calculate the 
parity violating asymmetry A in the ti system, additional renormalized amplitudes are 
needed. In terms of the tree- level amplitude, M , and the next-to-leading order SUSY 
QCD corrections, 5M, the renormalized amplitudes at the one- loop level can be written 
as M = Mo + 5M. Denote the momenta of the initial and the final state particles as 
1iij ) i)'im{p^) ~~ > ti(P2)ij(Pi), and the Dirac four-spinor as = u{j>i) (vi = v (pi)) for particle 
(anti-particle) i. Then, M = ig^T^T^Ji ■ J 2 /s, where Jf = v^p^^u^p^) and J 2 = 
u(p2) r y tJ 'v{pi); s is the invariant mass of the ti pair; g s and are the gauge coupling and 
the generator of the group SU(3) C , respectively. 

To calculate the parity violating asymmetry induced by the SUSY QCD effects, we fol- 



low the method presented in Ref . [0 , in which the asymmetry was calculated numerically 
using the helicity amplitude method. To obtain the renormalized scattering amplitudes, 
we adopt the dimensional regularization scheme to regulate the ultraviolet divergences and 



the on-mass-shell renormalization scheme [13| to define the input parameters. The SUSY 



QCD corrections to the scattering amplitudes arise from the vertex diagram, the gluon 
self-energy and the box diagrams, as well as the crossed-box diagrams. The renormalized 
amplitudes can be written as 

5M = 5M vl + 5M v2 + 5M S + 5M DB + 5M CB , (4) 

where SM vl and 5M v2 are vertex corrections, 5M S is the self-energy correction, and 5M 
and SM CB are the contributions from the box diagrams and crossed-box diagrams, respec- 
tively. The results for these separate contributions are, 

= ig 2 s {T^ m )u{p 2 ) • J x + Fffr + frff + fth + Ft • h 

+{FxVi + hFt 1 + rt'h + Ft 1 ■ Ji)7b)v(pi)/S, (5) 

5M^ = tgl{T^T[Jv{p,){Fffi 2 + rt-J2)u{p A )/s, (6) 
5M S = F S M , (7) 



6M DB = ^ 2 -(T^rf m )[F 1 ^(M 2 n 4 ^ 1 -M 27 5^375^i 
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+F 2 DB ^(-U 2 U 4 V 3 'J^V 1 - U^^UiVzVi - M 2 75M 4 ^37m75^1 + M27/.75M4^375^l) 
+ jpD fj,u (u^^UiVzTluVl + M27m75M4^37^75^i) 

+F® B (u 2 ^u i v ?J v 1 - 

J r F<P BtM (-U 2 r Y5U4V 3 J fM V 1 - U 2 ^^ b U A V^Vi - U 2 U 4 V^^ 5 Vi + U^^UaV^Vi) 

'(u 2 JfiU4y 37^75^1 + U2j,xj5U4V 3 J u Vi)}, (8) 
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6M CB = ^ s 2 -(T^r z c m )[Ff B (M 2M 3^i-M 2 7 5 ^475^) 

+F2 Bfl (-U 2 U 3 V i 'J fl V 1 - t^Ug^i - M275«3^47/.75^1 + ^27^75^475^1 ) 
+F CBp,u [u 2lfiU ^ ili/Vl + M 2 7 At 7 5 M 3 t;47 Jy 75fi) 

+F 4 C '' B (M275M3^1 - U 2 U 3 V^ 5 Vl ) 

+ ^? 5 -72275^3^47^^1 - M27m75«3^4^1 ~ ^3^47^75^1 + «27m U 3^475^i) 
+ ^ BMly (w 2 7 At W3^47^75^1 + M27^75M3^47^^l)], (9) 

where, fii = Ji^j^, etc., and the explicit expressions of the form factors , F i v \ Fq, 
, , and b i are given m the Appendix. 



Numerical Results and Conclusion 



In this section, we give our numerical results for a 176 GeV top quark jT4j. To avoid 



numerical instabilities and to take into account the fact that the decay products of the 
produced top quarks at large scattering angles are better distinguishable from the back- 
ground, we impose a cut on the transverse momentum (pt) and the rapidity (y) of the top 
quark and anti-quark: 

p T > 20 GeV and \y\ < 2.5. (10) 

We note that the parity violating asymmetry A would be independent of the parton dis- 
tribution functions (PDFs) if there were no kinematic cuts imposed in the calculation. 
Therefore, A should not be sensitive to the choice of the PDFs. In this paper we use the 
MRSA' PDFs |L5| and evaluate both the strong coupling a s and the PDFs at the scale 
Q = V~s = M tt -.Q 

As discussed in the previous section, the SUSY parameters relevant to our study are 
m tii m t2' ®t ( or m i L i m t R i m LR), m b R ' m qL an d m §- To simplify our discussion, we assume 
rriq L R = wig = m^ L , so that there are only four SUSY parameters to be considered, , 
mj 2 , 6i and m g . (The SU{2)l gauge symmetry requires that = m^.) 

The mSUGRA models predict radiative breaking of the electroweak gauge symmetry 
induced by the large top quark mass. Consequently, it is possible to have large splitting 
in the masses of the left-stop and the right-stop, while the masses of all the other (left- or 



right-) squarks are about the same ||17|| . For the MSSM models with scenarios motivated 
by current data @], a light ti is likely to be the right-stop (£r), with a mass at the order 
of mw\ the other squarks are heavier than t\. Since heavy superparticles decouple in loop 
1 Using the CTEQ4M PDFs 111] gives a similar result. 
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contributions, we expect that a lighter t\ would induce a larger asymmetry. Because the 
parity-violating effects from the SUSY QCD interactions arise from the mass difference 
between t\ and fa, it is obvious from Eq. ([[]) that the largest parity violating effect occurs 
when &i is ±7r/2 for m t - < rn^. When di = ±vr/4, the parity asymmetry should be zero. 
This is evident from the results shown in the Appendix, which indicate that the amplitudes 
that contribute to A are all proportional to Zi = =Fcos(26y. 

In either the mSUGRA or the MSSM models, the gluinos are usually as heavy as the 
light squarks, on the order of a few hundred GeV. However, Farrar has argued that light 
gluinos are still a possibility. If gluinos are light, then a heavy top quark can decay into a 
stop and a light gluino for < (m*— m§) such that the branching ratio of t — > hW + could 
show a large difference from that (~ 100%) predicted by the SM. The CDF collaboration 
has measured the branching ratio of t — > bW + to be 0.87ioioifui EH- At the la level, 
this implies that a 50 (90) GeV t\ requires the mass of the gluino to be larger than about 
120 (80) GeV for 6~ t = ±tt/2. However, at the 2a level (i.e. 95% c.L), there is no useful 
limit on the mass of the gluino.0 

To represent different classes of SUSY models in which the parity-violating asymmetry 
induced by the SUSY QCD interactions can be large, we show in Table 1 four represen- 
tative sets of models. They are labeled by the set of parameters (m^ , mi , 9i) , which are 
equal to (50,1033,-1.38), (90,1033,-1.38), (50,558,-1.25) and (90,558,-1.25), respec- 
tively. (All the masses are in units of GeV.) Based upon Eq. (2), one can also label these 
models by (m iL ,m iR ,m LR ), which are (1000,90,1100), (1000,118,1100), (500,40,520) and 
(500, 88, 520), respectively, for f3 = tt/4. 

It is interesting to note that for all the models listed in Table 1, the asymmetry A 
is negative (i.e. a^ < ol) for m g < 200 GeV, and its magnitude can be as large as 
3% for models with light t\. The maximal \A\ occurs when m g is about equal to (m t — 
77if ) because of the mass threshold enhancement. For m g > 200 GeV, the asymmetry 
A becomes positive, with a few percent in magnitude, and monotonically decreases as 
rrig increases. Comparing these results with those induced by the SUSY EW and SUSY 
Yukawa corrections H, it is clear that SUSY QCD interactions can generate a relatively 
larger parity-violating asymmetry. 

The differential asymmetry A{M t i) also exhibits an interesting behaviour as a function 
of the ti invariant mass M#. This is illustrated in Table 2 for the first SUSY model in 
Table 1 ((mj^m^,^-) = (50, 1033, -1.38)). As shown, |^4(M t j)| increases as M t i increases 
for m g < 200 GeV, which is similar to the effects from the SUSY EW and SUSY Yukawa 



2 A recent analysis on the possible experimental signature for a light gluino from top quark decay can 
be found in Ref. |20l. 
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Table V. 

Parity violating asymmetry A in pp — > ti + X, as a function of rrig, for four sets of SUSY 
models labeled by (m^m^, 



nig (GeV) 


(50, 1033, -1.38) 


(90, 1033, -1.38) 


(50, 558, -1.25) 


(90, 558, -1.25) 


2 


-1.10% 


-1.28% 


-0.98% 


-1.13% 


50 


-1.53% 


-1.77% 


-1.40% 


-1.59% 


100 


-2.34% 


-2.23% 


-2.21% 


-2.23% 


120 


-2.86% 


-1.99% 


-2.89% 


-1.99% 


135 


-3.16% 


-1.79% 


-3.43% 


-1.82% 


150 


-2.58% 


-1.50% 


-2.80% 


-1.53% 


175 


-1.18% 


-0.44% 


-1.30% 


-0.56% 


200 


0.99% 


0.97% 


0.82% 


0.77% 


225 


1.60% 


1.41% 


1.40% 


1.19% 


250 


1.53% 


1.34% 


1.35% 


1.17% 


275 


1.27% 


1.16% 


1.16% 


1.02% 


300 


1.04% 


0.94% 


0.95% 


0.83% 



contributions Q. When M t i is larger than 500 GeV, |^4(M t f)| can become very large at 
the expense of the signal event rate. For the model with rrig equal to 200 GeV, A{M t i) is 
positive for M t i less than ~ 500 GeV; it reaches its maximal value when M t i is about twice 
the gluino mass. This is again due to the mass threshold enhancement. 

We note that our conclusion for the magnitude of the parity-violating asymmetry is 
different from that given in Ref. |PH|.P| There the box (and crossed-box) diagram con- 
tributions were not included and some of the formulae (Eq. (8) and the form factor A 
in the Appendix) for calculating the vertex diagram contributions appear to contain mis- 
prints. From our calculation, we find that the box diagram contributions are in general 
small compared to the vertex diagram contributions, and are sizable only for heavy gluinos 
(rrig > 200 GeV). Therefore, the fact that we have included these corrections, while Ref. 
| f2~lf did not, does not constitute the main difference between the two calculations. To 
resolve this discrepancy, we compare our results with those in the literature for the SUSY 



QCD radiative correction to the production rate of qq — * tt at the Tevatron p2| . After cor- 
recting the relative sign between the box and crossed-box terms and the apparent misprints 
in the form factors i^f 3 and F^ 2 B Vi in Ref. |22| ,0 our results agree with those presented in 
[p2| , providing us with confidence in our results. However, our numerical results for the ti 
production rate do not agree with those given in Ref. . For the case considered in |H 



3 Without the cuts in (|Io|), the values of A for the first model in Table 1 are —1.0%, —2.65%, and 
-0.94% for rrig = 2, 120, 200 GeV, respectively. 



These apparent problems in Ref. [^2| were also pointed out in Ref. [^l) . 
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Table 2: 

The differential asymmetry A{M t i) and cross section da/dM t t (in unit of fb/GeV) as a 
function of M t j for the first SUSY model in Table 1 with various m g values. 



Mff (GeV) 

it V / 


m- g = 2 GeV 


m g = 120 GeV 


m-g = 200 GeV 


358 


-0.22% 


23.0 


-0.73% 


16.3 


44% 


19.4 


368 


-0.42% 


36.6 


-1.33% 


26.2 


0.95% 


31.4 


378 


-0.53% 


40.2 


-1.63% 


29.0 


1.39% 


34.9 


388 


-0.67% 


38.7 


-2.00% 


28.1 


2.12% 


34.2 


398 


-0.74% 


37.2 


-2.17% 


27.2 


3.33% 


34.4 


408 


-0.82% 


35.0 


-2.38% 


25.7 


3.67% 


32.4 


425 


-0.93% 


30.9 


-2.64% 


22.8 


2.35% 


27.0 


450 


-1.13% 


24.1 


-3.07% 


18.0 


1.05% 


20.0 


475 


-1.31% 


18.3 


-3.40% 


13.8 


0.11% 


14.7 


500 


-1.47% 


13.7 


-3.68% 


10.4 


-0.62% 


10.8 


525 


-1.61% 


10.3 


-3.81% 


7.9 


-1.21% 


8.0 


550 


-1.76% 


7.7 


-4.11% 


5.9 


-1.71% 


5.9 


575 


-1.93% 


5.7 


-4.34% 


4.4 


-2.08% 


4.3 


1000 


-3.26% 


0.032 


-5.82% 


0.026 


-4.66% 


0.024 



Table 3: 

The SUSY QCD corrections (Act) to the qq — > ti production rates at the Tevatron with 
s/s = 2TeV as a function of nig, for the first SUSY model in Table 1. 



m g (GeV) 


2 


50 


100 


120 


135 


150 


175 


200 


225 


250 


275 


300 


Aa (pb) 


1.17 


0.26 


-0.04 


-0.18 


-0.87 


-0.49 


-0.02 


0.33 


0.30 


0.24 


0.19 


0.16 



rrig = 200 GeV and mi = trig = 75 GeV, we obtain a 39%, in contrast to 33%, correction 
in the total cross section without cuts. Including cuts in (]i~0f) only slightly increases the 
correction to 40%. 

For completeness, in Table 3 we show the SUSY QCD corrections Act to the qq — > ti 
production rates at the Tevatron with ^/s = 2TeV, as a function of rrig, for the first 
SUSY model in Table 1. The SM tree level cross section, with the cuts in (10), is 4.06 
pb. At the upgraded Tevatron, it is expected to measure the cross section of ti pair 
production to ~ 10% (6%) with a 2 ftT 1 (10 fb" 1 ) integrated luminosity ||. Therefore, 
it is a challenging task to detect the difference in rates from SUSY QCD contributions 
predicted by these models unless gluinos are very light (of the order of 1 GeV). Nevertheless, 
the parity violating effect induced by the SUSY interactions could be observable with a 
large integrated luminosity (> 2 ftT 1 ) at the Tevatron. 
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Up to now, we have only considered the one loop SUSY QCD effects on the parity 
violating asymmetry A in ti pair production. Amusingly, the parity-violating asymmetry 
induced by the SUSY QCD interactions can also occur at the Born level. If gluinos are 
very light, of the order of 1 GeV, this asymmetry can be generated by the tree level process 
gg — > ti. Unfortunately, its production rate is smaller than the gg — > ti rate, which is only 
about one tenth of the qq — > ti rate at the Tevatron. Hence, it cannot be measured at 
the Tevatron. However, at the CERN Large Hadron Collider (LHC), the production rate 
of gg — > ti is large enough to allow the measurement of the parity-violating asymmetry 
induced by the SUSY QCD interactions. The asymmetry in the production rates of tjj 
and tni, generated by the gg fusion process alone, can reach about 10% for M# larger than 
about 500 GeV. We shall present its details and include the effect from the gg and qq fusion 
processes in a future publication |23 . 



This work is supported in part by the National Natural Science Foundation of China, 
and by the U.S. NSF grant PHY-9507683. 



Appendix 

We give here the form factors for the matrix elements appearing in Eqs. (8)-(12). 
They are written in terms of the conventional one-, two-, three- and four-point scalar loop 



integrals defined in Ref. [p4 



K X>1 = ^( 2m g Y i) Cfl (-P2,k,m- ti ,mg,mg) 



i=l,2 67T 
OL 

i=l,2 LZ7Y 

^ = £ ^[C 20 ^ - {m 2 t + ml)C Q Xi - 2m t m g C Y i \(-p 2 , k, m io m~ g , nig) 

i=l,2 ° 7F 

+ E TT-iBiXi - 2mtm s B' Yi + 2m 2 t B' 1 X l }(m 2 t , m h m~ ti ) 

i=l,2 671 
i=l,2 ° n 



Y ^{-2X i )Ci"(-p 2 ,k,mf l ,ms,m s ) 

i=l,2 ° n 
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^Z^C 20 + (m 2 - ml)C ){-p 2 , k,m~ u ,mg,m~ g ) 

i=l,2 871 " 

+ E 7T~ZiBi(m 2 , rrig, mi) 

i=l,2 37F 

53 -^ ± Z i m t C tx (-p 2 ,k,m i mg,rrig) 

i=l,2 87r 



£ ^(-2Z i )C7^(-p 2 ,fc,m^,mg,mg) 

i=l,2 S7F 

+ E ^[((P2-PirC'72-C'^)](-p 2 ^,m 5 ,m t -,m t -) 

i=l,2 iZ7F 

^(C 20 - m 2 gC )(p4, -k,mq,m~ g ,mg) 
2a 

47T 

+^[-(^2-^)^/2 - CH(p 4 , 
3a 



l r^Z-B 1 (m 2 q ,mg,mg) 



-2)C^ u (p A , -k, rrig, m~ g , m § ) 



({B 21 + Bi + ^~ (m 2 (B + 1) - 25 22 )/A; 2 )(A; 2 ,m § ,m § ) 
-(S 21 + £?! + i - m 2 B' + 2£? 22 ) (0, m 5 , m~ g )j 
+ E ^ ((25 22 (A; 2 , m 9 ~, m 9 ~) - A (m g -))/A; 2 - 25 22 (0, m,-, m 9 ~) 
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E ^((^ + m t) x i + 2mrngYdD (-p2,P4,P3,m i .,mg,m ij ,7 

i=l,2 47F 

E j 1 (mXi + mgY i )D^(-p 2 ,P4,P3,m it ,mg,mg,rrig) 

i=l,2 47r 

E T £x i^'^(-P2,P4,P3, mi., rrig, rrig, rrig) 

i=l,2 47F 

E ~7~Zi{m t — m~)Do(—p2,p4,P3, mi., rrig, rrig, rrig) 

i=l,2 47r 

53 -j-Z i m t D iX {-p 2 ,p4, p 3 , m~ u ,m~ g , m q , rrig) 

i=l,2 47F 



i=l,2 47F 



In the above: 

k = pi +P2 = P3 +P4, s = k 2 , 

Xi = aj + b\ = 1, YJ = - b 2 = T sia(20 f ), Z< = 2aA = q= cos(20 t -). 
where ai, 61, a 2 , b 2 are given by 

ai = -^=(cos0 f -sin0f), 61 = --^=(cos0 t ~ + sin0 t ~), 

02 = Wcos0 t ~ + sin^) = 61, b 2 = /=(cos0£ — sin^-) = — ai. 

v2 v2 



Also, 



, _ dB (p 2 ,mi,m 2 ) , _ dB^p 2 , m u m 2 ) , _ dB 22 (p 2 ,mi,m 2 ) 

g p2 , *1~ g p2 ~ , *22~ - g p2 

c 20 = 9lUf c^ - \. 
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